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Abstract. The high-frequency integrated spectra of 
Compact Steep Spectrum (CSS) sources show breaks with 
a moderate spectral steepening well fitted by continu- 
ous injection synchrotron spectra. In lobc-dominated CSS 
sources the radiative ages deduced by the synchrotron the- 
ory are in the range of up to 10^ years, if equipartition 
magnetic fields arc assumed. These radiative ages are well 
correlated with the source size indicating that the CSS 
sources are young. In order to maintain the frustration 
scenario, in which the sources' lifetimes are « 10^ years, 
their equipartition magnetic field would be systematically 
decreased by a factor ^20. To complete the sample used 
in this work, we conducted observations at 230 GHz with 
the IRAM 30-m telescope of those sources which did not 
have such high-frequency observations up to now. 

Key words: Radio continuum: galaxies; radiation mech- 
anisms: non-thermal 



1. Introduction 

Flux density limited radio source samples are known to 
contain a large portion of compact sources with high fre- 
quency spectral index a > 0.5 (S^ ^ v^"^) and angular 
sizes below 2". The corresponding projected linear sizes 
are typically < 15kpc^. Fanti et al. (1990) have shown 
that the majority of these CSS sources cannot be larger 
sources foreshortened by projection effects, which means 
that their radio emission originates on sub-galactic scales. 

It is obvious that one should investigate the connec- 
tion between CSS and larger radio sources. Two scenarios 
have been proposed that would naturally explain the ob- 
served small sizes. First, they could reflect an early stage 

Send offprint requests to: M. Murgia 

^ We assume Ho = 100 km s~^ Mpc~^ and qo = 0.5 



in the evolution of radio sources. This is the youth sce- 
nario (Phillips & Mutel 1982; Carvahlo 1985). The second 
possibility is that the unusual conditions in the interstel- 
lar medium of their host galaxies, such as a higher density 
and/or the presence of turbulence, inhibit the radio source 
from growing to larger sizes. This is the frustration sce- 
nario (van Breugel et al. 1984). 

Most CSS sources exhibit double-lobed structures such 
as seen in classical radio galaxies. This feature is com- 
mon to both quasars and galaxies. Quite a few of these 
are symmetric, which gave rise to the terms CSO (com- 
pact symmetric objects, with sizes < 0.5 kpc) and MSO 
(medium-size symmetric objects, with sizes > 0.5 kpc). 
These sources are considered as scaled-down versions of 
larger-sized double radio sources (Fanti et al. 1995). A 
minority of CSS sources is made up by sources with com- 
plex or highly asymmetric structures. These are mostly 
quasars, with most of their luminosity provided by the 
jet. The distortions may be caused by jet bending, which 
is further amplified by strong projection. 
A plausible organization scheme for these different mor- 
phological classes is to identify CSO and MSO with pro- 
genitors of large doubles ('baby radio galaxies'), while the 
asymmetric CSS sources could represent frustrated radio 
sources. 

2. Sample selection 

The sample we choose for this study consists of: 

1) the sample selected from the 3CR and the Peacock & 
Wall (1981; PW) catalogues (Fanti et al. 1995), and 

2) nine sources extracted from the B3VLA sample (Vi- 
gotti et al. 1989). 

1) The 3CR/PW sample (38 sources) was established 
on the basis of the well-known criteria: linear size < 15 
kpc, spectral index a > 0.5 and a limit in radio power 
P > 1026-8 W/Hz at 178 MHz for the SCR, or P > lO^^ 
W/Hz at 2.7 GHz for the PW. High resolution structure 
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information is available for these sources from the liter- 
ature. In the framework of this project wc observed at 
230 GHz the 17 sources for which the flux densities at 
this frequency were not available from literature. 

2) The second complete sample of CSS sources was se- 
lected from the "strong" section of the B3 VLA catalogue 
('S'408MHz > 0.8 Jy) according to the following definitions: 
'^limiHz 0.5 up to 10 GHz, and a size of 15 kpc. 
Sources with 0:458^^2 < 0.5 have also been included if 
their spectrum is steeper than a = 0.5 a.t v > 1.4 GHz. 
These are Giga-Hertz-pcaked spectrum (GPS) candidate 
sources. This sample contains 83 objects. The sample has 
been recently observed with the VLA at 5 and 8.4 GHz, 
with resolutions of ^ 0'.'2. Additional observations with 
VLBI arrays are planned. Only about 60% of these sources 
have at the moment either a photometric or a spectro- 
scopic redshift, therefore we limited the list to the objects 
which could be observed at 230 GHz to the latter. Due to 
scheduling constraints we only observed 9 sources of the 
B3VLA sample. 

3. Observations and data reduction 

The observations were done with the IRAM 30-m tele- 
scope near Granada, Spain between 10th Feb. and 15th 
Feb., 1998. The receiver we used was the 37-channel 
bolometer described by Kreysa et al. (1998), with feeds 
arranged in a hexagon around the central feed, and beam 
separations of 20". The sensitivity of each channel was 
70mJy/s. The beam size was 10'.'5 as derived from point- 
ing scans. All our sources were point-like to this beam. 
Opacity measurements were made via skydips, from which 
we derived zenith opacities of 0.1-0.8. Despite the high 
atmospheric attenuation of the signal, the observations 
(generally done at high source elevations) could be per- 
formed satisfactorily, owing to the relatively quiet atmo- 
sphere present during our observing run (resulting in a 
relatively low sky noise). The calibration factor to con- 
vert the observed counts into flux densities was determined 
by observing the source K3-50A, which has a total flux 
density of 10 Jy and a peak flux density of 7.6 Jy for the 
beam size of our measurements. The measures were finally 
corrected for atmospheric attenuation and the elevation- 
dependent gain of the telescope. The pointing was fre- 
quently checked by cross-scanning the sources 0923-1-392, 
1144-1-402, 1308-1-326, and 1418-1-546. The pointing accu- 
racy was found to be ^ 3". The same sources were used 
to adjust the focus at regular intervals. 

The sources were observed in an ON-OFF mode in 
which the wobbling secondary mirror of the telescope 
moved at a frequency of 2 Hz between the source po- 
sition and positions located at ± 60" in azimuth. Each 
source was observed between 3 and 5 times, with 10 ON- 
OFF pairs each. The differences of the flux densities be- 
tween the main and the reference beams were averaged by 
weighting inversely proportional to the square of the r.m.s. 



noise of each measurement. The resulting integration time 
spent on each source was thus typically 50 minutes, and 
the final r.m.s. noise of each measurement was between 2 
and 5 mJy. In addition, the calibration may be uncertain 
by 15% owing to day-to-day variations in the system gain. 

The data were recorded with all receivers, with the 
central one pointing at the source. As the angular extents 
of the target sources are in the range of a few arcseconds, 
the outer channels contain only sky emission. Hence, they 
were used to subtract a mean sky level from the central 
channel. Using this procedure, we assumed that the scale 
size of the fluctuations of the sky emission is large com- 
pared to the total coverage of the bolometer array (60"), 
which corresponds to 1.5 m at a distance of 5 km. Together 
with the wobbling rate, this ensures efficient rejection of 
sky noise, since the time scale of atmospheric fluctuations 
is much larger (see e.g. Altenhoff et al. 1987). 

In Tab. 1 we list the measured flux densities and 1-a 
errors, inclusive of noise and calibration uncertainty, at 
230 GHz. 



Table 1. Flux densities at 230 GHz 



Name 


Pointing centre 


(B1950) 


S 


230 








a 






S 


[mJy] 




3C-PW 






















0404-1-76 


04 


04 


00 


13 


76 


48 


52.50 


115 


± 18 


0428+20 


04 


28 


06 


86 


20 


31 


09.13 


62 


± 11 


0740+38 3C186 


07 


40 


56 


77 


38 


00 


30.97 


5 


± 


4 


0758+14 3C190 


07 


58 


45 


05 


14 


23 


04.40 


8 


± 


3 


1005+07 3C237 


10 


05 


22 


04 


07 


44 


58.56 


14 


± 


5 


1019+22 3C241 


10 


19 


09 


38 


22 


14 


39.63 


2 


± 


3 


1153+31 


11 


53 


44 


08 


31 


44 


16.00 


9 


± 


3 


1203+64 3C268.3 


12 


03 


54 


08 


64 


30 


18.40 


1 


± 


2 


1225+36 


12 


25 


30 


77 


36 


51 


47.00 





± 


3 


1358+62 


13 


58 


58 


36 


62 


25 


06.71 


29 


± 


6 


1416+06 3C298 


14 


16 


38 


80 


06 


42 


21.30 


8 


± 


4 


1447+77 3C305.1 


14 


47 


49 


35 


77 


08 


46.65 


2 


± 


3 


1517+20 3C318 


15 


17 


50 


63 


20 


26 


52.95 


5 


± 


3 


1634+62 3C343 


16 


34 


01 


06 


62 


51 


41.80 


12 


± 


4 


1637+62 3C343.1 


16 


37 


55 


29 


62 


40 


34.24 


7 


± 


4 


1819+39 


18 


19 


42 


33 


39 


41 


15.00 


5 


± 


4 


1829+29 


18 


29 


17 


93 


29 


04 


58.29 


9 


± 


4 


B3VLA 






















0809+404 


08 


09 


31 


62 


40 


28 


02.76 


11 


± 


3 


0810+460B 


08 


10 


58 


61 


46 


05 


48.12 


2 


± 


2 


1025+390B 


10 


25 


49 


34 


38 


59 


57.55 


46 


± 


8 


1128+455 


11 


28 


56 


40 


45 


31 


24.57 


4 


± 


3 


1159+395 


11 


59 


16 


35 


39 


35 


52.89 


5 


± 


3 


1225+442 


12 


25 


15 


70 


44 


17 


17.23 


1 


± 


3 


1233+418 


12 


33 


10 


88 


41 


53 


38.00 


11 


± 


4 


1242+410 


12 


42 


26 


39 


41 


04 


29.97 


83 


± 13 


1350+432 


13 


50 


24 


04 


43 


14 


09.40 


4 


± 


3 
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4. Results cind discussion 

For each soiircc we have compiled flux densities at different 
frequencies from the literature, mostly from Kiihr et al. 
(1979) and from the CATS database (Verkhodanov et al. 
1997). Our new measures at 230 GHz have been added to 
the compilation. All flux densities have been brought to 
the BGPW scale (Baars et al. 1977). Examples of source 
spectra are shown in Fig. 1. For the fit algorithm described 
in Sect. 4.2 we have assumed flux densities less than '3a 
to be upper limits only. 

Most of the sources show significant departure from the 
classical power law which describes a zero age transparent 
synchrotron spectrum from a relativistic electron popula- 
tion with power law energy distribution. The deviations 
from the power law are of the following type: a) a low- 
frequency turnover (the most conspicuous deviation); b) a 
steepening at high frequencies. High-frequency flattening, 
if any, is quite rare. The above deviations are interpreted 
as due to synchrotron self-absorption and to particle en- 
ergy losses, respectively. In order to describe them, we 
have fitted the compiled flux densities with a synchrotron 
aged spectrum Sagod(i^) (described in the next section), 
modified by low-frequency absorption, as follows (Pachol- 
czyk 1970): 

S{U) OC {u/u.r+f'il - e-(-/''i)"'°+'') . Saged(i-) (1) 

where i>i is the frequency at which the optical depth is 
equal to 1. In case of an homogeneous synchrotron self- 
absorbed source /3 = 2.5, while a is the not aged spectral 
index in the transparent frequency range. 



4.1. The synchrotron aged spectrum model 

We assume that the radio source evohition is described by 
a continuous injection model, where the sources are contin- 
uously replenished by a constant flow of fresh relativistic 
particles with a power law energy distribution, with expo- 
nent 6. It is well known that, under these assumptions, the 
radio spectrum has a standard shape (Kardashev 1962), 
with spectral index ajnj = {6 — l)/2 below a critical fre- 
quency fbr and ah — ainj -1-0.5 above vy,,-. If there is no ex- 
pansion and the magnetic field is constant, the frequency 
i/br (in GHz) depends on the elapsed time since the source 
formation, Tgyn (in Myrs), the intensity of the magnetic 
field B (in /xG) and the magnetic field equivalent to the 
microwave background Bomb = 3.25(1 + zf' (in [iG) as: 

" '''°B^ + BgMB K(1 + .)]V2 (2) 

The whole spectral shape cannot be described by an an- 
alytic equation, the two behaviours described being only 
the asymptotical ones, and has to be computed numeri- 
cally. This model is referred as continuous injection (CI). 
Fitting the spectral data to the numerically computed CI 
spectrum, one obtains the break frequency i4>r, from which 
the source age is obtained if the magnetic field is known. 
This simple model does not consider expansion effects, 
which may be important if the source is young. So, the 
simplest variant of the original model is the one in which 
the radiating particles loose energy through expansion and 
the magnetic field changes according to flux conservation 
(Kardashev 1962, CIE). An alternative possibility to be 
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considered, in an expanding source, is that the magnetic 
field changes less rapidly than in the flux conserving as- 
sumption because of a continuous magnetic flux input as- 
sociated to the fresh particles injection. We set B oc 
where for m = 2 wc get the flux conserving expansion. 
This is referred as CIm model. In our special case we as- 
sume m = 1, consistent with the models applied by Bald- 
win (1982) or Begclman (1996). Although the theoretical 
background to these models is in the Kardashev paper, 
we have decided to present in the Appendix the detailed 
development. The spectral shapes of the CIE and CIm 
models have been computed again numerically. The break 
frequency for the CIE model is sixteen times higher than 
in the CI case, for equal elapsed time and final magnetic 
field intensity (Kardashev 1962). In the CIm model, in- 
stead, one finds that the break frequency is (2.5m — 1)^ 
larger than in the CI model. The asymptotic behaviours 
at frequencies lower and larger than the break frequency 
are the same as in the CI model. However the steepening 
occurs over a broader frequency interval. In order to better 
emphasize the differences between these three models it is 
useful to consider, together with the usual flux- frequency 
representation, the point-to-point spectral index defined 



CI models 

TTTI| 1 I I lllll| 1 I I lllll{ 1 I I lllll{ 1 I I lllll{ 1 I I lllll{ 1 I I lllll{ 1 I I lllll{ 



as 



dlog Saged(i^/'^br) 
rflog U/Ubr 



(3) 



Both representations are shown in Fig. 2. While in the 
flux-frequency plane the differences are hardly visible, 
they can be much better traced in the point-to-point spec- 
tral index behaviour. Note that the displayed figures only 
show the theoretical difFcrcnccs. In practice the observed 
spectra permit only fits in the flux-frequency plane. 

4-2. The spectral fits 

Spectral fits to the spectra have been made with the CI, 
CIE, and CIm models. They allow us to determine the 
non-aged spectral index ainj and the break frequency Vbr, 
that, together with the normalization, are the three free 
parameters characterizing all the models. The best fits are, 
surprisingly, obtained with the CI model. The other mod- 
els have steepenings which are too gradual for the majority 
of the spectra. An example is given in Fig. 3. The reduced 
Xred of the models including adiabatic expansion and a 
variable magnetic field is always greater (typically twice) 
than that of the CI model (Fig. 4). The CI fits appear 
quite good, even in cases of high x^ed values, which, at vi- 
sual inspection, appear more due to an under-estimation 
of the flux density errors than to a poor fit of the spectral 
model on the data. 

The majority of the spectra show a clear break fre- 
quency, with a change of slope Aa 0.5. We stress that 
there is no evidence for spectral steepening with Aa sig- 
nificantly larger than 0.5. Only a few sources are fitted by 
simple power laws. In these cases u^r could be either very 




Fig. 2. Tlie three continuous injection models described in the 
text. The solid line corresponds to the CI, the dashed line to the 
CIE and the dot-dashed line to the CIm. On the top the flux 
density (arbitrary units) is plotted as a function of the ratio 
v/vbr (flux- frequency plane), on the bottom the point-to-point 
spectral index is shown as a function of v/iAjt- The injection 
spectral index is the same for all three models ainj = 0.625. 



high (^ 100 GHz) or very low (^ few 100 MHz). In these 
sources wc have preferred the low frequency choice, since 
for the high frequency one would have implied abnormally 
high values for ajnj as compared to the other sources. 

Wc also note that only a very small fraction of the CSS 
sources, if any, shows some evidence of fiux density excess 
at high frequency, as it would be caused by a flat spectrum 
core or by thermal dust emission. Perhaps the only case 
is 3C138, where the core, known from VLBI observations, 
shows up in the integrated spectrum at 230 GHz only. 

It may appear surprising that the flts with the assumed 
spectral model are so good. In fact the sources of the 
sample consist of several components, as lobes and jets 
and hot spots, where physical conditions can differ from 
one another and therefore also break frequencies may be 
different. It is likely that the spectrum is dominated by 
the brighter component(s). In addition, one could imag- 
ine that some confusion might have occurred between gen- 
uine spectral steepening due to energy losses and the low- 
frequency turnover due to absorption processes. We feel 
that this is a minor problem, but it is difficult to quantify 
it (see, however, next section). 

The break frequencies range from a few hundred MHz 
to tens of GHz. At low frequencies, as said above, the limit 
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Table 2. Physical Parameters 



Source 


redshift 


LS 




Qinj 




Tsyn 




type 








[kpc] 


[GHz] 




[10"'^G] 


[10' yrs] 






3C-PW 




















0127+23 


(3C43) 


1.46 


9.40 


>100 


0.76 


5.0 


0.3 


55.0 


b 


0134+32 


(3C48) 


0.37 


2.10 


1.8 


0.50 


2.0 


11.3 


0.60 


b 


0138+13 


(3C49) 


0.62 


3.60 


2.6 


0.57 


7.0 


1.3 


4.39 


b 


0221+67 


(3C67) 


0.31 


6.80 


2.2 


0.56 


0.7 


51.0 


0.22 


a 


0223+34 




1.00 


3.80 


24.1 


0.38 


9.4 


0.3 


24.20 


b 


0316+16 




1.00 


1.20 


12.0 


0.81 


9.8 


0.3 


5.74 


b 


0345+33 


(3C93.1) 


0.24 


1.20 


1.6 


0.58 


1.0 


35.8 


0.05 


? 


0404+76 




0.60 


0.53 


10.7 


0.38 


4.9 


1.1 


0.76 


a 


0428+20 




0.22 


0.45 


7.2 


0.38 


5.8 


1.2 


0.59 


b 


0429+41 


(3C119) 


1.02 


0.75 


9.5 


0.49 


8.0 


0.5 


2.36 


b 


0518+16 


(3C138) 


0.76 


2.90 


5.3 


0.47 


1.0 


16.7 


0.28 


b 


0538+49 


(3C147) 


0.55 


2.40 


1.6 


0.44 


2.5 


8.3 


0.47 


b 


0740+38 


(3C186) 


1.06 


8.20 


0.3 


0.75 


0.7 


113.2 


0.12 


a 


0758+14 


(3C190) 


1.20 


14.10 


30.3 


0.79 


0.7 


10.6 


2.15 


b 


1005+07 


(3C237) 


0.88 


4.50 


2.4 


0.59 


1.2 


18.1 


0.40 


a 


1019+22 


(3C241) 


1.62 


2.80 


1.1 


0.77 


2.2 


9.3 


0.49 


a 


1153+31 




1.56 


3.20 


6.9 


0.69 


1.7 


5.4 


0.95 


a 


1203+64 


(3C268.3) 


0.37 


3.90 


4.4 


0.66 


0.7 


35.4 


0.18 


a 


1225+36 




1.98 


0.17 


3.9 


0.64 


18.0 


0.2 


1.41 


b 


1250+56 


(3C277.1) 


0.32 


4.40 


0.7 


0.39 


0.4 


204.3 


0.04 


a 


1323+32 




0.37 


0.18 


6.9 


0.38 


5.1 


1.4 


0.20 


a 


1328+30 


(3C286) 


0.85 


14.20 


7.3 


0.38 


13.0 


0.3 


77.82 


b 


1328+25 


(3C287) 


1.06 


0.40 


10.7 


0.47 


3.5 


1.7 


0.39 


b 


1358+62 




0.43 


0.16 


71.9 


0.70 


9.0 


0.2 


1.39 


b 


1416+06 


(3C298) 


1.44 


9.10 


<0.1 


0.50 


1.6 


>50 


<0.3 


a 


1443+77 


(3C303.1) 


0.27 


5.00 


0.8 


0.64 


0.6 


110.2 


0.07 


a 


1447+77 


(3C305.1) 


1.13 


9.00 


0.8 


0.62 


0.4 


151.3 


0.10 


a 


1458+71 


(3C309.1) 


0.90 


7.80 


109.4 


0.65 


7.0 


0.6 


21.0 


b 


1517+20 


(3C318) 


0.75 


7.80 


4.8 


0.69 


1.5 


9.6 


1.32 


b 


1607+26 




0.47 


0.18 


6.3 


0.71 


10.0 


0.5 


0.55 


a 


1634+62 


(3C343) 


0.99 


0.70 


3.8 


0.65 


2.0 


6.5 


0.17 


? 


1637+62 


(3C343.1) 


0.75 


1.30 


1.7 


0.62 


1.5 


16.0 


0.13 


a 


1819+39 




0.40 


3.10 


3.3 


0.73 


6.0 


1.6 


3.13 


b 


1829+29 




0.60 


9.30 


12.1 


0.63 


6.0 


0.8 


19.12 


b 


2248+71 


(3C454.1) 


1.84 


6.60 


0.4 


0.69 


1.5 


^26.5 


<0.40 


a 


2249+18 


(3C454) 


1.76 


2.10 


11.9 


0.72 


5.0 


0.8 


4.29 


b 


2252+12 


(3C455) 


0.54 


13.70 


1.7 


0.62 


0.2 


348.2 


0.06 


a 


2342+82 




0.74 


0.66 


10-100 


0.79 


4.5 


^1.3 


0.82 


a 


B3VLA 




















0809+404 


0.55 


4.80 


4.8 


0.53 


1.1 


16.1 


0.48 


? 


0810+460B 


0.33 


2.80 


13.9 


0.94 


1.1 


10.3 


0.44 


a 


1025+390B 


0.361 


5.10 


>100 


0.65 


0.8 


6.0 


1.40 


a? 


1128+455 


0.40 


1.90 


1.8 


0.53 


1.5 


17.3 


0.18 


? 


1159+395 


2.37 


0.35 


2.8 


0.38 


4.7 


1.6 


0.35 


? 


1225+442 


0.22 


0.50 


2.5 


0.67 


0.8 


40.3 


0.02 


? 


1233+418 


0.25 


2.70 


3.5 


0.51 


0.6 


52.7 


0.08 


a? 


1242+410 


0.811 


0.40 


4.7 


0.38 


2.6 


4.2 


0.16 


? 


1350+432 


2.149 


8.00 


0.2 


0.84 


1.2 


55.4 


0.23 


? 
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Fig. 3. Three diflterent fit results for tlie source 2252+12. Tlie 
solid line corresponds to a pure CI model, the dashed line rep- 
resents the CIE model, the dot-dashed line stands for the CIm 
model. Note the different qualities of the fits, expressed by the 
various x^cd values! 

is set by confusion with the effects of absorption processes. 
The injection spectral index ainj ranges from 0.35 to 0.8, 
with (ainj) — 0.63. The typical errors of the break fre- 
quencies and the injection spectral indices as given by the 
fit algorithm are up to 40% and 0.05, respectively. The 
results of the CI fits are compiled in Tab. 2. 

4-3. Radiative ages and the nature of CSS sources 

In order to determine radiative ages, from Eq. (2) and 
variants of the other models, the magnetic field B has to 
be known. We stress that the age depends rather strongly 
on B, which is somewhat uncertain. We take the equipar- 
tition value Bcq. Our assumption is motivated by the 
fact that Beq accounts rather well for the low-frequency 
turnover in terms of synchrotron self-absorption. We are 
aware that this is not a proof for equipartition. Other au- 
thors (e.g. Bicknell et al. 1997) prefer instead thermal ab- 
sorption. In any case, Beq represents a poor statistical 
upper limit to B, in the sense that, were it larger by a 
factor of four, the low frequency turnovers would be sys- 
tematically higher than observed by « 30%. 

Using the value of Seq, we obtain, from Eq. (2), ra- 
diative ages Tgyn ranging from 10^ to 10^ years. Since 
the intrinsic magnetic fields of the CSS sources in our 
sample strongly overweight the magnetic field equivalent 
to the cosmic microwave background (-Bcmb)7 the latter 



Fig. 4. Xicd histograms for the CI, CIE, and CIm models. On 
the bottom right the histogram of the ratio between the CIm 
and the CI Xrcd is also shown. 

can be neglected in Eq. (2). Therefore, if the source mag- 
netic field deviates by a factor / from the field determined 
for equipartition Bcq, the radiative ages will change by 
j-3/2_ These radiative ages do not necessarily represent 
the source ages, but rather the radiative ages of the dom- 
inant source component(s). Only when the lobes, which 
have accumulated the electrons produced over the source 
lifetime, dominate the source spectrum, the radiative age 
Tgyn is likely to represent the age of the source. If, in- 
stead, the spectrum is dominated by a jet or by hot spots, 
the radiative age likely represents the permanence time of 
the electrons in that component and is expected to be less 
(perhaps much less) than the source age. In addition, dom- 
inant jets or hot spots might have their break frequency 
up-shifted by Doppler effects. 

The existing structure information on our sample, 
mostly from MERLIN and VLBI observations, allows us 
to split the sources in two groups: those in which the over- 
all spectrum is dominated by lobes (classified "type a" in 
Tab. 2); those in which the spectrum is dominated by a 
bright jet or hot spot (classified "type b" in Tab. 2). For 
the B3 VLA sample the available information does not yet 
allow such a morphological sub-division. 

The sources of class b have radiative ages systemati- 
cally lower than those of class a, as we expected. Further- 
more, while the synchrotron age is well correlated with 
the source size for class a, it seems that there is no corre- 
lation at all between the linear size and the radiative age 




Fig. 5. Linear size as a function of tlie synclirotron age for type 
a (filled dots) and type b (open dots) sources. The B3VLA 
sources have been excluded. The horizontal dashed lines rep- 
resent the selection limits of the linear size distribution of the 
sources in our samples while the diagonal lines reflect constant 
values of I^HiE.. 



for class b sources (see Fig. 5). We have further computed 
for each source the expansion velocity Woxp ~ LS/2rsyn, 
where LS is the source largest dimension. The two classes 
show very different distributions: 
Class a (wcxp/c) = 0.34 ± 0.07 cr„„^p = 0.28 
Class b (woxp/c) = 11 ± 4.4 cr,,^^^ = 20 
Provided the assumed magnetic field is reasonably correct, 
the above values for the radiative ages indicate that the 
CSS sources are young. We further note that the ages, and 
corresponding expansion velocities, are not far (somewhat 
larger) from the recent results on expansion of CSOs by 
Owsianik et al. (1998) and Owsianik & Conway (1998). 
As the radiative ages are strongly dependent on the as- 
sumed magnetic field, a field only a factor two lower than 
assumed would be required for a better agreement. 
In order to maintain the frustration scenario, in which 
the sources' lifetimes are « 10'' years, their equipartition 
magnetic field should be decreased by a factor ^ 20. 
One could think that the correlation between linear sizes 
and the radiative ages shown in Fig. 5 could be a par- 
tial consequence of the equipartition assumption. In fact, 
Bcq oc LS^^^^ implies that LS oc rj/n • This is not the case 
for the following reasons: 1) no correlation between linear 
sizes and radiative ages is obvious for type b sources. 
2) In particular, the break frequencies seem to be corre- 
lated with the linear sizes for type a sources (Fig. 6). 



Fig. 6. Linear size as a function of the break frequency (in 
the source rest frame) for type a (filled dots) and type b (open 
dots) sources. The B3VLA sources have been excluded. The 
horizontal dashed lines represent the selection limits of the 
linear size distribution of the sources in our samples. 

These means, at least for class a sources, that the break 
frequency is an effective clock indicating the source age. 
The correlation seen in Fig. 5 is not an artifact completely 
introduced by the assumption of equipartition in Eq. (2). 

5. Conclusions 

The high-frequency integrated spectrum of CSS sources 
shows a break with a moderate spectral steepening (Aa « 
0.5), well fitted by the continuous injection model spec- 
trum with constant magnetic field. Spectral fits with syn- 
chrotron models with decreasing magnetic field are def- 
initely poorer. There is no evidence of sharp cut-offs as 
would be expected if the supply of relativistic electrons 
had stopped more than « 10* years ago. 

In lobe dominated CSS sources the radiative ages Tgyn 
are likely to represent the source ages. Assuming equipar- 
tition magnetic field strengths the source ages are in the 
range of up to 10'^ years, in agreement with the recent 
results on expansion velocities of small size CSOs by 
Owsianik et al. (1998) and Owsianik & Conway (1998) 
derived from VLBI observations. Therefore, assuming the 
CSS sources to be young, the magnetic field has been de- 
duced to be equal to the equipartition value within a factor 
« 2. 

In jet or hot spot dominated CSS sources, the radiative 
life times are much shorter and very likely represent the 
permanence time of the electrons in those components. 
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Appendix 

A. Synchrotron radiation, adiabatic expansion and 
continuous injection in the case of a time- 
Vciriable magnetic field 

We assume that adiabatic expansion and injection of par- 
ticles start at time to- The source is continuously replen- 
ished by a constant flow of particles with a power law en- 
ergy spectrum with spectral index S intensity B decreases 
according to 



B 



and the characteristic size R of the regions occupied by the 
particles (which should not be confused with the source 
size LS) grows as 



R = Ro 



The continuity equation for this model is 

dN d , n 5 , t 



dt 



= ai 



(Al) 



N{e, t) is the diff'erential energy spectrum of the relativis- 
tic particles at time t. The first term on the right hand side 
describes the effects of the synchrotron losses, ai = CiBj_, 
ci is constant (Pacholczyk 1970). The second term ac- 
counts for adiabatic expansion (02 = r/t). The term 03 
describes the injection rate and is considered to be con- 
stant for simplicity. 
The solution of Eq. (Al) is 



N{e,t) 



,\r(S-l) 



5-2 



cIt' 



(A2) 



where t\ov, is MAX[fo, Tmin] and Tmin is implicitly given by 



- 1 



ciBH 



2m+r-l 



The energy e^j.{T) is 
* ( ^ 2m + 



m+r— 1 



1) 



e^j.(T) is the break energy of the particle population in- 
jected at the time t { to < r < t). The solution (A2) 
represents the sum of all the populations injected form 
T = to to T = t. The integral (A2) can be written as 



N{e,t) =a3e-^-'^t 



2m + r-l 



f 



(1-2) — 



dz 



(A3) 
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where ziow is equal to MAX[0, 1 — e/ebr 
2m + r - 1 



ebr 



and 



CiBlt 



In the energy spectrum two break energies are always 
present, ebr and ebrmin) the break energy of the first pop- 
ulation injected att = to. These regions can be identified 
in the energy spectrum: 

1) e < ebrmin • 

In this region the spectrum is a power law: 



r{S - 1) + 1 



1 - 



r((5-l) + l' 



2) ebrmin ~ e < Cbr- 

The spectrum starts to gradually deviate from the low 
energy power law at e « ebrmin and undergoes a steep- 
ening at e PS Cbr- 



N{e,t) 



ebr 



2m -|- r 

ebr 



(1-z) 



3) e > Cbr- 

The spectrum is again a power law but with a spectral 
index steeper by 1. 

-5-1 



036 



ciBi(5-i) 



Since the separation between these two energies increases 
with time the shape of the spectrum also changes with 
time. But an asymptotic stationary shape still exists when 
t/to — > 00. In fact, in this case ebrmin 0, and at every 
energy e, z\o^ = 0. In this way the integral 



N{e,t) = ase-^-H 



ebr 



/" 

Jo 



2m + r-l 



{1-zf-f + l 



z'-'dz 



does not depend on the ratio t/to anymore. 

The integral can be solved only by numerical means, but 

analytic asymptotic limits can still be found. 



N{e,t) 
N{e,t)'. 



r{S - 1) -h 1 



e < 



e > 



2m + r-l 

ciB\t 
2m -t- r — 1 
ciB\t 



(A4) 



i.e. in the spectrum only the Cbr break energy is present. 
The asymptotic solutions are two power laws. The low- 
energy power law spectral index is 5, the high-energy 



power law spectral index is 5-1-1. Moreover, the normaliza- 
tion for e ^ ebr is time-independent, i.e. a perfect balance 
is reached between the number of particles that leave this 
region because the synchrotron and expansion losses and 
the number of particles constantly injected into the source. 



A.l. CI model 

The CI model reproduces the simple situation in which 
there is no adiabatic expansion and the magnetic field 
strength stays constant. This corresponds to r = and 
TO = in (A3). In this case the integral can be directly 
solved resulting in: 



N{e,t) = 
N{e,t) = 
where 

ebr = 



036' 



-<5-l 



ciBl{5-l) 



1 - (1 - e/ebr) 



5-1 



03 e 



-5-1 



ciBl{5-l) 



e < ebr 



e > ebr 



CiBi(t-to) 

The shape of the spectrum does not depend on time (Kar- 
dashev 1962). 

A. 2. CIE model 

The CIE model reproduces the situation in which the vol- 
ume containing the particles is adiabatically expanding 
at a constant rate and the magnetic field is frozen in the 
plasma. Since R oc t, the conservation of the magnetic 
flux requires that B oc fr"^ . This corresponds to r = 1 and 
m = 2 in (A3). The shape of the spectrum changes with 
time during its rise. However an asymptotic stationary so- 
lution exists when f:^ to: 



N{e,t) 



ebr = 



ciBl Jo 



(1-z) — 



1 



" (5-2 I 

z dz 



ciB\t 



In this asymptotic limit the CIE spectrum shows a break 
energy four times greater than the break energy of the CI 
model. Using (A4) one finds that 



N{e,t) 

N{e,t) 



036-^ 



-5-1 



ciBl{5-l) 



e < 



e » 



cxB\t 

4 

ciB\t 



this solution is the same found by Kardashev (1962). The 
asymptotic solutions are two power laws with the same 
dependence on the spectral index as in the case of the CI 
model. However, the energy range necessary to complete 
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the transition between the two asymptotic power laws is 
wider. For t ^ to the break energy of the CIE model is 



distributed and that the time scale for their continuous 

re-isotropization is much shorter than the radiative time- 



greater than the break energy ebrmin of a CI model with scale, ci{Bj_) = gCiS (Jaffe & Perola 1974). 
the same final magnetic field. This result appears a bit 
surprising at first glance. In fact, due to the expansion References 
losses and the stronger magnetic field, the break energy of 
the first population of particle injected at to in the CIE 
model is lower than that of the corresponding population 
of the CI model. However, both the stronger mean mag- 
netic field and the expansion losses will also decrease the 
energy of the first populations in a way that, at time t, 
they do not contribute appreciably to the integrated spec- 
tra. 



A. 3. CIm model 

The CIm model reproduces the situation in which the vol- 
ume containing the particles is adiabatically expanding 
and there is always equipartition between magnetic field 
and particle energy density. The equipartition condition 
u oc {u is the energy density of the relativistic parti- 
cles) implies that R oc f™/^. This corresponds to r = m/2 
in (A3). An asymptotic stationary solution exists when 
t~> to: 



N{e,t)^ 



ebr 



ciBl 

2.5m - 1 
ciBlt 



(1-z) 



1 



z'-^dz 



In this asymptotic limit the CIm spectrum shows a break 
energy (2.5m — 1) times greater than the break energy of 
the CI model. Using (A4) one finds that 



N{e,t) 
N{e,t) 



036 



■1) + 1 

-(5-1 



ciBl{S-l) 



€ > 



2.5m- 1 
ciBlt 

2.5m- 1 
ciBH 
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For the CIm the CIE comments are also valid. 



A. 4-. Emission spectrum 

The emission spectrum is given by the convolution of the 
emission spectrum of the single electron F [v / {C2B sinOe^)) 
(with C2 constant, Pacholczyk 1970) with the energy dis- 
tribution N{e, t): 

Sa.ed(./^.,ainj)0C^^^ F • 

N{x ■ ebr,i) dx dO 

where z^r = C2-Bebr^ and ainj = {5 — l)/2. It is as- 
sumed that the pitch angles 6 between the electron ve- 
locity and the magnetic field direction are isotropically 
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